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The US Navy’s GEOSAT Follow-On spacecraft (GFO) primary mission objective is 
to map the oceans using a radar altimeter. Satellite laser ranging data, especially in 
combination with altimeter crossover data, offer the only means of determining high-quality 
precise orbits. Two tuned gravity models, PGS7727 and PGS7777b, were created at NASA 
GSFC for GFO that reduce the predicted radial orbit hrough degree 70 to 13.7 and 10.0 mm. 
A macromodel was developed to model the nonconservative forces and the SLR spacecraft 
measurement offset was adjusted to remove a mean bias. Using these improved models, 
satellite-ranging data, altimeter crossover data, and Doppler data are used to compute both 
daily medium precision orbits with a latency of less than 24 hours. Final precise orbits 
are also computed using these tracking data and exported with a latency of three to four 
weeks to NOAA for use on the GFO Geophysical Data Records (GDR’s). The estimated 
orbit precision of the daily orbits is between 10 and 20 cm, whereas the precise orbits have 
a precision of 5 cm. 


Nomenclature 


Subscripts 

l spherical harmonic degree 

m spherical harmonic order 

Conventions 

Cim,Si m Normalized spherical harmonic coefficients 
Pim Normalized associated Legendre function 

c,i Drag coefficient 

GM Gravitational constant, m 3 /s 1 2 
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&2 Love number of degree two 
kz Love number of degree three 
r radius 

R e Reference radius for gravity model 
U Gravity potential 

CDDIS Crustal Dynamics Data Information System 
DGFI Deutsches Geodatisches Forschungsinstitut 
GFO Geosat-Follow-On 

GOTOO.2 Goddard Ocean Tide Model 2000, Version 2 

GPS Global Positioning System 

IGDR Intermediate Geophysical Data Record 

ILRS International Laser Ranging Service 

IRI International Reference Ionosphere 

LRA Laser Retroreflector Array 

MOE Medium precision Orbit Ephemerides 

NCEP National Center for Environmental Prediction 

opr once-per-revolution 

POD Precision orbit determination 

POE Precise Orbit Ephemerides 

SLR Satellite Laser Ranging 

SWH Significant Wave Height 

T/P TOPEX/Poseidon 

WVR Water Vapour Radiometer 

Symbols 

(3 Specular reflectivity 
S Diffuse reflectivity 
A longitude 

</> latitude 


I. Introduction 


Table 1. GFO Orbit Summary 


I. A. Background 

The GEOSAT Follow-On (GFO) spacecraft was launched on February 10, 1998 with the objective of providing 
continuous ocean observations along the GEOSAT exact repeat ground track for both real-time and near- 
real-time measurements of relative ocean heights, and absolute heights for large-scale ocean modelling. The 
inclination and ground-track repeat period (17 days) complement the data collected by other missions, such 
as TOPEX/Poseidon, ERS-1, ERS-2, ENVISAT and Jason-1 (see Table 1). 

GFO carries a single-frequency (13.5 Ghz) altimeter, a dual- 
frequency water vapour radiometer, a dual-frequency Doppler 
(TRANET-style) beacon for operational tracking, and a laser retro 
reflector array (LRA) for precision orbit determination (see Fig. 1). 
The satellite also carried GPS dual-frequency receivers, however the 
GPS system on GFO only supplied limited data, and could not be 
used for precision orbit determination (POD). Hence, both the oper- 
ational and precise orbits have been determined using a combination 
of satellite laser ranging (SLR) and Doppler tracking in combination with the use of altimeter crossovers. 
The spacecraft was declared operational on November 29, 2000, and has now returned nearly six years of 
altimeter data over 120 repeat cycles. GFO has occupied the GEOSAT Exact Repeat Orbit and supplied 
useful data longer than the original GEOSAT mission which provided data in this orbit from 1986 to 1989. 1,2 

The GFO spacecraft was constructed by Bail Aerospace and launched on a Taurus launch vehicle from 
Vandenburg Air Force Base on February 10, 1998 for the prime customer, the U.S. Navy. The National 
Oceanic and Atmospheric Administration has the responsibility for distribution of the altimeter data. The 
Planetary Geodynamics Branch at the NASA Goddard Space Flight Center produces the operational and 
precise orbits for GFO. The precise orbits are supplied to NO A A for placement on the Geophysical Data 


Altitude 
Eccentricity 
Inclination 
Repeat Period 


800 km 
0.008 
108° 

244 revs in 17 days 
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Records (GDR’s) which are released to the scientific community. 

I.B. Orbit Determination Challenge 

For altimeter satellites, the prime observation is the range mea- 
surement from the satellite to the ocean surface. It follows that 
the ability to discriminate changes in the height (or topogra- 
phy) of the ocean depends on the on knowledge of the satellite 
orbit, and that the satellite orbit must be known as accurately 
as possible. The precision of the orbit depends on the quality 
of the tracking data, the fidelity of the force and measurement 
models, and the choice of parameters estimated in an orbit de- 
termination solution. It is the knowledge of the radial compo- 
nent of the orbit that is the most critical for satellite altimetry. 

Over the 13 years of its mission, T/P achieved radial orbit pre- 
cision of 2-2.5 cm. 3,4 For the successor mission, Jason-1, orbits 
with a radial precision of 1-cm have been achieved. 4,5 Both 
T/P and Jason-1 benefited directly from the near-continuous 
tracking supplied by DORIS (Doppler Orbitography and Radio 
Positioning Integrated by Satellite) or GPS (Global Position- 
ing System). Indirectly, both missions benefited from years of effort to improve the gravity models and the 
nonconservative force modelling, both major components in the orbit error budget. The situation with GFO 
is more nuanced. The lack of continuous tracking makes the achievement of radial precision comparable to 
Jason or T/P more difficult. The challenge then becomes how we are to use the available tracking (SLR, 
Doppler, altimeter crossovers) to achieve the highest quality orbits. In the absence of the GPS data from 
GFO, we also had to demonstrate that the SLR data could be used to produce operational orbits with a 
latency of less than 24 hrs. In this paper we describe the quality of the satellite tracking data, the improve- 
ments to the force modelling that have been implemented, and the principal factors that affect the GFO 
radial orbit precision. 



II. Data 


II. A. Satellite Laser Ranging 

The international network of satellite laser ranging stations operates under the aegis of the International 
Laser Ranging Service (ILRS). 6 This network is shown in Fig. 2. The network is dominated by stations 
in the Northern Hemisphere, most notably in Europe. The satellite laser ranging data can have a precision 
of a few mm, especially for the best stations in the network. In reality, of the 40 or so stations that are in 
the network, only a subset provide data on a routine basis. In Fig. 3, we show the number of passes of 
satellite laser ranging data acquired by GFO between January 1, 2005 through March 2006. A satellite pass 
is defined as a single passage of GFO over a tracking station for which data were acquired. Tracking data 
are only acquired when the station is staffed, and when the weather permits optical tracking. An additional 
consideration is the tracking priority assigned to each satellite target by the ILRS. Generally, the satellites 
lowest in altitude receive the highest priority. The list of tracking priorities is reevaluated on a regular basis, 
based on scientific needs or changes in the operational status of the user satellites. In this priority scheme, 
GFO as an altimeter satellite at relatively low altitude (800 km) receives a high ranking. 

From Fig. 3, the Australian station, Yarragadee (YARA), is the prime contributor of SLR data for 
GFO and supplies 15% of the data. The remainder of the stations in the NASA network (Monument Peak, 
California [MNPE], Greenbelt, Maryland [GRF1], and McDonald Observatory, Texas [MCDOl]) supply 11% 
of the passes for GFO. The rest of the international network supplies 74% of the SLR data for GFO, with the 
largest contributors being Zimmerwald (ZIMM, Bern, Switzerland), Graz (GRAZ, Austria), Herstmonceux 
(RGO, East Sussex, U.K.), Wettzell (WETT, Germany), and Riyadh (RIYA, Saudi Arabia). The network 
provides, on average, between 12 to 14 passes of SLR data per day. 
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Figure 3. Number of SLR tracking passes of GFO ordered by station, from January 2005 through March 2006 
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II.B. Doppler Data 


The Doppler data from GFO are obtained from three stations: Guam, Point Mugu, California, and Prospect 
Harbor, Maine. The Doppler data are in the style of data from the OPNET stations that tracked GEOSAT, 
and consist of one-way Doppler data transmitted from GFO, and received at the Earth tracking stations. 1,2,7 
The data have a noise of approximately 2 cm/s. Although this is quite noisy by modern standards, the 
Doppler data still help to tie down the orbit when the SLR tracking is sparse. In addition, the Doppler data 
are needed to obtain new orbit solutions following orbit trim maneuvers. 

II. C. Altimeter Data 

In order to supplement the SLR and Doppler tracking, we use the altimeter data supplied by GFO as an 
additional data type. The data are used in the form of altimeter crossovers. 8,9 Altimeter crossover data 
are formed by differencing altimeter ranges from two intersecting passes and interpolating to the point of 
interesection. 10 The GFO altimeter crossovers provide dense spatial coverage over the ocean areas. The 
data are edited in regions with high sea surface variability (greater than 20 cm), and in shallow seas (less 
than 500 m depth). In addition, a maximum residual edit criterion of 20 cm is also applied. The altimeter 
range corrections, as applied in the IGDR (Intermediate Geophysical Data Record), are listed in Table 2. 
We use a T/P derived tide model, GOTOO.2, for the ocean tide altimeter range correction. The GOTOO.2 
tide model is based on 286 cycles of altimeter data from TOPEX/Poseidon. The GOTOO.2 tide model is an 
update of the GOT99.2 tide model, which was based on 232 cycles of TOPEX/Poseidon altimetry. 11 The 
dry troposphere correction is derived from the NOAA Global Forecast System, based on the Global Data 
Assimilation System. 13 For the wet troposphere correction we use first the GFO water vapour radiometer 
(WVR) correction. If the IGDR contains a null field for the GFO WVR, we apply the NCEP model derived 
value. We tested the IRI2000 ionosphere model, 14 but did not discern a significant improvement with respect 
to using the IRI95 15 model for GFO. 


Table 2. Altimeter Range Modelling for the GFO IGDR 


Ocean tide 

GOTOO.2 

Earth tide 

Cartwright and Eden (updated) 12 

Dry troposphere 

NCEP 13 

Wet troposphere 

GFO WVR or NCEP 13 

Ionosphere 

IRI95 15 

Inverse barometer 

/(dry troposphere) 16 

EM bias 

3.8% SWH 


III. The GFO Orbit Determination System 


III. A. Overview 

The orbit determination system we have designed for GFO imports tracking data and ancillary data from 
a variety of sources, and delivers three products: the medium precision orbits (MOE’s), the precise orbits 
(POE’s), and the ephemeris predictions for the satellite laser ranging stations. The satellite laser ranging 
data are delivered at least once daily from the tracking stations to the ILRS data centers, at the NASA 
GSFC Crustal Dynamics Data Information System (CDDIS) in Greenbelt, Maryland and at the Deutsches 
Geodatisches Forschungsinstitut (DGFI) in Munich, Germany. 17, 18 The data centers exchange their recently- 
delivered data daily to ensure that their holdings are equalized so that users can obtain reliably obtain data 
should one of the data centers be temporarily inaccessible. The GFO Doppler data are delivered electronically 
from the Naval Space Operations Center (NAVSOC) at Point Mugu, to NASA GSFC Monday through Friday, 
but not on weekends. The IGDR altimetry data are obtained daily from NOAA. On a typical day, the SLR 
and Doppler data are imported by early afternoon local time ( 16:00 to 17:00 hrs UT). This means we 
use SLR and Doppler passes that are obtained through the day of the arc, even up to mid-afternoon UT 
time. Each MOE arc spans five days ending on the current day, but due to latency of processing at NOAA, 
altimeter crossovers are included in the MOE’s only for the first three days of the MOE arc. 
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In addition to the tracking data, we must import other ancillary information from external sources. We 
obtain daily updates on Earth orientation from the International Earth Rotation Service (IERS) at the Paris 
Observatory. 19 Solar flux and geomagnetic data are downloaded daily from NOAA National Geophysical 
Data Center. 20 The F10.7 solar flux and geomagnetic indices are used as inputs to atmospheric density 
models which are required to model atmospheric drag at the GFO altitude. 

The MOE orbit is exported daily by 21:00 hrs 
to 23:00 hrs UT to data users. The precise orbit 
(POE) is computed with a latency of three to four 
weeks. The MOE orbit is computed daily using a 
sliding window of five days of data. The POE arcs 
are six days in length and overlap by only one day. 

The POE’s benefit from stabilized values for the 
Earth orientation, and use altimeter data to form 
crossovers over the entire data arc. The daily MOE 
arc is used to create an ephemeris prediction which 
is supplied to the laser stations to ensure continued 
tracking for GFO. 

Orbit trim maneuvers disrupt normal satellite 
operations and require special attention. The orbit 
trim maneuvers maintain the satellite to within ± 

1 km of the nominal GEOSAT ground track. How- 
ever they invalidate the predicted ephemerides sent Figure 4. Near-Real-Time Orbit Determination System for 
to the laser tracking stations during the previous 
days. In 2005, there were 11 orbit trim maneuvers. 

The SLR stations must be notified of the impending maneuver in order to prevent fruitless attempts to 
track the satellite. NAVSOC supplies post-maneuver state vector predictions for GFO. We propagate this 
state vector and distribute these predicted ephemerides to the stations to ensure continued tracking of the 
satellite. 

Spacecraft anomalies also perturb normal orbit system operation. If the anomalies cause the spacecraft to 
deviate from its nominal nadir attitude, the SLR stations must also be notified since the laser retroreflector 
array might not oriented properly to permit SLR tracking. 

We use the NASA GSFC GEODYN orbit determination and orbit parameter estimation program to 
compute the GFO orbits. This orbit determination program is a batch, least-squares filter that processes 
numerous types of satellite tracking data. 21 GEODYN is also used to compute the orbits of T/P and Jason- 
1, 4 and processed the tracking data that was incorporated into the gravity models developed to support 
GFO. 22,23 A schematic of the GFO orbit determination system is depicted in Fig. 4 The analysis of the 
GFO tracking data requires detailed force and measurement modelling which we will summarize in the next 
sections. 

III.B. Force Modelling 

We must model as accurately as possible all the forces that affect the spacecraft trajectory. For altimeter 
satellites, it is the gravity field induced error, and the nonconservative force model error that are the largest 
contributors to the radial orbit error budget. 3 

III.B, 1. Gravity Field 

We model the gravity field in spherical harmonics using normalized coefficients ( Ci m , Si m ) using the equa- 
tion 24 



m (sin <fl) [Cim cos(mA) + S lm sin(mA)] (1) 

r 1=2 771= o ' ' 

where GM is the universal constant of gravitation times the mass of the Earth, l is the degree, m is the 
order, Pi m are the fully normalized associated Legendre polynomials, R e is the reference radius of the Earth, 
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4> is the latitude, and A is the longitude. By definition, the degree one terms are zero, since we choose the 
origin of the coordinate system to be at the center of mass of the Earth. 
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Figure 5. Radial orbit error vs. spherical harmonic order for GFO from the gravity field error covariances through 
degree and order 70. 


We currently use the PGS7777b gravity model for GFO orbit 
production. The model, developed at NASA GSFC, was based on 
CHAMP and other satellite data, and included data based 54 arcs 
of GFO tracking data in 2000 to 2001, including SLR, Doppler, and 
altimeter crossovers. Effectively, this model was tuned for GFO and 
used as a base the GPS tracking data and precision accelerome- 
try data from CHAMP. 23 The projected radial orbit error was re- 
duced from 65.2 mm with JGM-2 25 to 49.8 mm with the JGM- 
3 gravity model, 26 and 26.2 mm with EGM96 27 to 10 mm with 
PGS7777b. The JGM-2 and JGM-3 gravity models were developed 
for the TOPEX/Poseidon mission, so that although they contained 
GEOSAT tracking data, they were not specifically tuned to the 
GEOSAT orbit. EGM96 included a complete reprocessing of the JGM-2 satellite tracking data, contri- 
butions from other satellites, and most importantly much improved altimetry and surface gravity data. 
EGM96 included GEOSAT Doppler and altimetry data from November 1986 through January 1987 but the 
the contribution of the GEOSAT data was limited by the short span of data, and the lack of direct ties 
between the TRANET and OPNET Doppler stations that tracked GEOSAT and the SLR data that in large 
part defined the low degree field and reference frame for EGM96. 

The PGS7727 model, developed at GSFC as a derivative of EGM96, included both TOPEX/GFO (dual- 
satellite) and GFO/GFO (single satellite) altimeter crossovers, including the SLR data to GFO. The addition 
of these data reduced the radial orbit error to 13.2 mm. 22 We show the gravity field error as a function of 
spherical harmonic order in Fig. 5. The JGM-3 and EGM96 gravity models had large errors at order one, and 
at the resonant and the near-resonant orders (e.g., m =14,15,28,29). GFO has heightened sensitivity to the 
gravity field at order one due to the m-daily perturbations. 24 It is clear that the GEOSAT Doppler tracking 
alone in these gravity solutions was insufficient to model adequately these gravity field terms for GFO. The 


Table 3. Gravity Field Radial Orbit 
Error for GFO 


Gravity model 

Orbit error 
(mm) 

JGM-2 (1993) 

65.2 

JGM-3 (1996) 

49.8 

EGM96 (1996) 

26.2 

PGS7727 (2001) 

13.2 

PGS7777b (2003) 

10.0 
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addition of GFO SLR and altimeter crossover data in PGS7727 and PGS7777b was highly effective in tuning 
the gravity field at those orders. 

The PGS7777b includes secular variations in C 20 , C 21 , and S 21 , and annual variations <$ 620 , SC 30 and 
( 5 ( 5,10 (see Table 4). The periodic variations in the low degree zonals are derived from an independent analysis 
of SLR and DORIS tracking data spanning two decades . 28 


Table 4. PGS7777b Time- Variable Gravity- 
Terms 


Term 

Value’ 


C 20 (Secular) 

1.258 x 

10 

-11 

/yr 

C 21 (Secular) 

-0.337 x 

10 n /yr 

£21 (Secular) 

1.606 x 

10 

-11 

/yr 

5(520 Annual Cosine 

3.4163 

X 

10 - 

-11 

< 5 C 2 o Annual Sine 

10.2115 

X 

10 

-11 

<5(530 Annual Cosine 

-0.3148 

X 

10- 

-11 

<5(530 Annual Sine 

0.1877 

X 

10- 

-11 

< 5 C' 4 o Annual Cosine 

-3.6574 

X 

10- 

-11 

<5(540 Annual Sine 

-4.1864 

X 

10- 

-11 


* The epoch of the PGS7777b solution is 1998.0. 


III.B.2. Macromodel 

The solar radiation pressure, planetary radiation pressure (due to the Earth’s albedo and thermal emission), 
and the atmospheric drag are modelled using a macromodel . 29 Analagous to TOPEX, the GFO spacecraft 
is modelled as a series of elemental flat plates, and the contributions of each plate are vectorially summed 
to obtain the total acceleration (See Fig. 6 ). 




Figure 6. Schematic of macromodel for GFO. 


We model GFO using 8 elemental plates representing the +X, -X, -Y, +Z, -Z spacecraft faces, the front 
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side of the solar array, and front and the back of the radar altimeter. Self-shadowing complicates the 
development of the macromodel. The +Y face of the spacecraft is obscured by the solar array and does 
not enter into the computations. Since the solar array tracks the Sun, the -Z surface are shadowed. Hence, 
in the macromodel, the -Z reflectivity parameters are set to zero. Thus, the -Z plate participates in the 
drag calculations, but not in the radiation pressure calculations. We list the macromodel areas, and the 
specular (/3) and diffuse (5) reflectivities in Table 5 which are used in the current generation of MOE and 
POE orbits. The GFO macromodel does not account perfectly for self-shadowing effects that might vary 
with beta prime (the angle of the Sun above or below the orbit plane). In addition, the macromodel does not 
include radiation interactions between surfaces. Also, unlike T/P, a detailed thermal model for POD was 
not developed. Notwithstanding these imperfections, the GFO macromodel is an improvement over using a 
simple cannonball for the nonconservative force modelling. 


Table 5. GFO Macromodel Parameters 


Plate 

Area 

(m 2 ) 

Reflectivities 

(A 

TX face 

0.320 

(0.36, 0.09) 

-X face 

0.736 

(0.44, 0.11) 

-Y face 

2.370 

(0.526, 0.132) 

+Z face 

2.450 

(0.555, 0.139) 

-Z face 

0.750 

(0.0, 0.0)* 

Solar array (front) 

3.987 

(0.144, 0.04) 

Radar altimeter (front) 

0.880 

(0.688, 0.172) 

Radar altimeter (back) 

0.880 

(0.688, 0.172) 


* The -Z plate surface area is reduced to account for shad- 
owing by the solar array. In addition, this plate partic- 
ipates in the macromodel drag calculations, but not the 
radiation pressure calculations. 


III.B.3. Other Force model effects 

The solid earth tides are modelled with ^2=0. 300 and £ 3 = 0 . 093, and special modelling for the free core nuta- 
tion. 27 The ocean tides use as a background the GOT99.2 ocean tide model, 1 1 derived from TOPEX/Poseidon 
altimetry. The PGS7777b gravity solution also included adjustments for the resonant tide terms (as described 
in Ref. 27) and these are overlain on the background tide model. We use the MSIS86 atmospheric density 
model. 30 The planetary radiation pressure due to the Earth’s shortwave and longwave flux is modelled as 
described in Ref. 31, where the total acceleration at each time step is computed as a double summation 
over both the contribution from each plate and the elements on the surface of the Earth in view of the 
satellite. The third body perturbations due to the Sun, Moon, and planets are modelled using the DE403 
set of planetary ephemerides. 32 

III.C. Measurement Modelling 

III.C.l. Station coordinates 

The SLR station coordinates are based on the ITRF2000 solution. 33 Coordinates for new SLR stations that 
were not originally part of the ITRF2000 solution (for example the new station at Mt Stromlo, Australia, 
rebuilt after the destruction of the old station there due to bushfires in January 2003) were obtained from 
ground survey ties, or independent adjustments using data to the satellites Lageos-1 and Lageos-2. The 
TOPEX POD team applied a number of corrections to the coordinates of some of the SLR stations in 
ITRF2000. The SLR station coordinates used for GFO are identical to those used in the production of the 
precise orbits for the TOPEX/Poseidon geophysical data records. NAVSOC supplied a priori coordinates 
for the Doppler stations which were adjusted to the SLR frame. Ocean loading corrections are computed 
using the GOTOO.2 ocean tide model and include the M2, S2, K2, Kl, 01, PI, Ql, Mf, Mm, and Ssa tidal 
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constituents. We include both the vertical, as well as the eastward and northward displacements due to ocean 
loading. The ocean loading correction diminishes as the distance of the station from the coast increases. 
The ocean loading corrections must be included in order to exploit the precise SLR data. For example, for 
the M2 constituent, the vertical displacement load tide might be 10 mm or less for the Mt. Stromlo SLR 
station, but 20 to 30 mm for the San Fernando SLR station. 

III.C.2. Measurement Offset Modelling 

It is necessary to define two vectors in the spacecraft 
coordinate system: (1) the location of the spacecraft 
center-of-mass (CoM), and (2) the location of the 
phase center of the laser retroreflector array. While 
the LRA is fixed, the CoM varies with propellant 
usage (see Fig. 7). Early in the mission, the ini- 
tial analyses of SLR data showed a large mean in 
the residuals. The SLR data were used to estimate 
the LRA offset. We cannot say whether the adjust- 
ment was a results of an error in the LRA location, 
or an error in the specification of the CoM in the 
spacecraft coordinate system. Nonetheless, the ap- 
plication of the new offset reduced the SLR residual 
mean from -2.5 cm to zero, over the test period. The 
offset adjusted by -6.3 cm in X, -1.1 cm in Y, and 
+11.1 cm in Z. 

III.C.3. Attitude Modelling 

The GFO spacecraft follows a yaw steering algo- 
rithm that maintains the solar array pointed at the 
Sun, while the radar altimeter is pointed at nadir. Unlike TOPEX/Poseidon, where quaternions are routinely 
available during off-nominal attitude regimes, for GFO we must rely completely on an analytical description 
of the spacecraft attitude. We obtained a limited amount of quaternion data for GFO early in the mission 
and compared the orientation angles for GFO computed from the analytical attitude model with those from 
the satellite telemetry. As depicted in Fig. 8, the differences are on average less than 0.3° in roll, pitch, and 
yaw. Thus, as long as GFO follows its prescribed attitude law, we will be able to orient the macromodel 
correctly in inertial space. In addition, assuming a maximum 0.3° error in attitude, we are assured that the 
error in the LRA measurement offset correction due to incorrect attitude knowledge will be less than 4 mm. 

IV. Results 



Figure 7. GFO laser retroreflector array geometry in the 
spacecraft coordinate system. 


IV. A. MOE orbits 

For the MOE orbits, we adjust the spacecraft state, a drag coefficient ( Cd ) per day, and empirical once- 
per-revolution (opr) accelerations along-track and cross-track to the orbit once per arc. Range-rate and 
troposphere biases are adjusted for the Doppler data. Pass-by-pass biases are adjusted for some of the 
non-core SLR stations of the ILRS. The Doppler data are weighted at 2 cm/s, whereas the SLR data are 
weighted at 10 cm. However, the data from some SLR stations are downweighted to between 60 cm and 1 
meter. The altimeter crossover data are weighted at 10 cm. All arcs with altimeter crossover data adjust 
an altimeter timing bias whose typical value is between 0.8 to 1.3 msecs. The adjustment of a timing bias 
on the altimeter data makes the time tag on the altimetry consistent with the time recorded at the SLR 
stations. 

Prior to January 20, 2004, the PGS7727 gravity model was used for MOE production. The PGS7777b 
gravity model has been used since January 20, 2004. The altimeter crossovers were introduced into the 
solutions for the MOE orbits starting on February 4, 2004. The crossovers were introduced to strengthen the 
solution for the MOE orbits, and compensate for diminished tracking schedules at the NASA SLR stations 
due to budget issues. 
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Figure 8. Comparison of GFO roll, pitch, and yaw from telemetry and from analytical model. 


We show the RMS of fit for the daily MOE arcs in Fig 9. The long-term average daily RMS of fit is 6.1 
cm for the SLR data, and 7.3 cm for the crossover data. However, there are large variations, and there is 
a correlation with the F10.7 solar flux with the RMS of fit degrading during high flux periods. The RMS 
of fit is also typically higher for arcs that start after maneuvers. The GFO orbit maneuvers almost always 
cause a disruption in the SLR tracking at some level, even if the predicted post-maneuver ephemerides are 
distributed in a timely fashion. In addition, for reasons of latency, for the first 2-3 days, the daily post- 
maneuver arcs do not contain any altimeter crossover data. Thus, the MOE arcs in the immediate aftermath 
of maneuvers have less dynamical strength than the ’normal’ five-day MOE arcs. 

Spacecraft anomalies will also affect the quality of the MOE orbits. Sometimes the cause of the degraded 
MOE orbit quality is not always obvious. For example, from November 4 to November 13, 2004, the RMS 
of fit to the SLR data in the daily arcs increased to between 13 and 25 cm. The precise orbits also show an 
increase in the SLR RMS of fit at that time, though only to 9 cm and only for the six-day arc beginning 
on November 7, 2004. It is likely there was a satellite or a data glitch on approximately November 10 or 
11, 2004. Because of the sliding window used for the MOE’s, if a satellite or a data anomaly is not properly 
identified, then the effects can be felt in the MOE’s for several days. 

Another way to assess orbit consistency is through the use of orbit overlap comparisons. We routinely 
compare the orbits from sequential MOE arcs, as long as there is no interruption due to a maneuver. Hence, 
the orbit overlaps usually contain four days of common data. We must be cautious in interpreting these 
statistics as they will be optimistic in their orbit quality assessment since they have so much data in common. 
We show these daily overlaps in Fig. 10 and summarize the full statistics in Table 6. The mean RMS overlaps 
are 5.11 cm radially, 11.89 cm cross-track, and 23.82 cm along-track. We omit the first arcs after maneuvers 
or spacecraft anomalies which would otherwise distort the results. We see once again, that maneuvers always 
cause a disruption in orbit consistency, and that it takes several days for the orbit statistics to stabilize. A 
semiannual signal is evident in the radial and cross-track orbit overlaps which might be related to errors in 
modelling of the atmospheric drag or the radiation pressure. 

IV. B. POE Results 

The POE orbits are computed with a latency of three to four weeks. The extra latency allows extra SLR 
and altimeter data to be imported to cover more of the data arc. In addition, the POE processing can 
take advantage of any independent analysis of the SLR or altimeter data that points to station, data, or 
spacecraft anomalies that were not diagnosed for the MOE production. Finally, the latency permits us to 
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Days Since January 1 , 2004 

Figure 9. SLR and crossover RMS of fit for daily MOE (medium precision orbit) arcs from February 4, 2004 through 
July 28, 2006. 





(a) Radial (b) Cross-track (c) Along-track 

Figure 10. RMS orbit overlaps for daily sequential MOE orbits, radial, along-track and cross-track to the orbit. The 
orbit overlaps are usually 4 days in length, except after orbit trim maneuvers. 
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Table 6. Orbit Overlap Statistics for Daily MOE arcs from Febru- 
ary 4, 2004 to July 28, 2006 


MOE Orbit Statistics (cm)* 



Mean 

Median 

Minimum 

Maximum 

Radial 

5.11 

2.49 

0.38 

479 

Cross-track 

11.89 

42.39 

0.43 

730 

Along-track 

23.82 

97.62 

1.78 

1521.6 


* The orbit statistics are computed over 836 daily MOE arcs. Arcs 
on the first day after a maneuver or a spacecraft anomaly are 
omitted from this summary. 


import stabilized and updated values for the Earth orientation and the solar flux. The set of estimable 
parameters is the same as for the MOE arcs, except for the empirical accelerations and drag coefficients 
(cds). 

Since altimeter data are available for the full six days of the data arc, more crossovers can be computed, 
especially since the number of potential crossovers is proportional to n 2 , where n is the number of orbital 
revolutions in a data arc. A typical POE arc has on average 2600 crossovers, whereas a typical MOE arc 
has only 675 crossovers over the first three days of the five-day arc. The extra data mean that a denser 
parameterization can be employed to accommodate force model error. Hence, for the POE’s we apply the 
same parameterization as for the TOPEX POE’s: we adjust opr’s along-track and cross-track once per day, 
and we adjust drag coefficients (c</ s) every eight hours. Between January 4, 2000 and June 28, 2006, we 
computed and exported 423 POE arcs. At first, we computed the GFO POE’s with the PGS7727 gravity 
model, however once the PGS7777b gravity model became available in early 2004, all the earlier POE arcs 
were recomputed with the newer gravity model. For the PGS7777b generation POE’s, the average SLR RMS 
of fit is 4.37 cm and the average crossover RMS of fit is 7.51 cm. We show the RMS of fit in Fig. 11 and 
summarize the RMS by year in Table 7. The RMS for all data types (SLR, crossovers, Doppler) peaks in 
2001 and 2002. The F10.7 radio flux from the Sun peaked between January 2000 and early 2002, so these 
peaks in the GFO RMS of fit are roughly coincident with the peak of solar cycle 23. The RMS of fit of a 
GFO POE arc may vary with the solar cycle for two reasons: (1) At the peak of the solar cycle, GFO is more 
susceptible to mismodelling from atmospheric drag; (2) For a single- frequency altimeter, the ionosphere error 
is more substantial near the solar cycle peak, than near the solar minimum. 

We also examine the orbit overlap statistics for the POE arcs (see Table 8), and compare these to the 
MOE orbit overlap statistics (see Table 6. From these statistics, we can guess that the radial orbit precision 
is 2.5 times the mean radial orbit overlap of 1.84 cm. 

In order to make an orbit accuracy assessment for GFO, we can estimate the sea surface variability 
using the GFO orbits, and then repeat the calculation after adjusting the GFO orbits using TOPEX/GFO 
crossovers. We illustrate this process in Fig. 12. We first computed the GFO sea-surface variability using the 
PGS7727 orbits and then adjusted the GFO orbits empirically relative to TOPEX/Poseidon. The resultant 
sea surface variability is 8.3 cm, compared to 9.5 cm before adjustment. The RSS difference between the sea 
surface variability maps is 4.6 cm, and can be taken as an estimate of the radial orbit error for GFO (relative 
to TOPEX). This analysis was done using the PGS7727 orbits early in the mission, so we have confidence 
that the GFO POE’s now computed with the PGS7777b gravity model have a radial precision of 5 cm. 

V. Conclusion 

The GFO mission was rescued by the on-board presence of the laser retroreflector, and the demonstration 
of near-real time POD using a combination satellite laser ranging data, Doppler and altimeter crossover data. 
Starting in January 2000, we have produced a 6.5 year span of POE orbits for the GFO GDR. With these 
orbits, altimeter analysis has shown that GFO is a Poseidon-class altimeter, and provides an important 
supplement to Jason- 1, TOPEX/Poseidon (T/P), ERS and ENVISAT. GFO altimeter data are being used 
to monitor inland lakes, and has been used for near-real-time monitoring of hurricanes and to detect the 
2004 Boxing Day tsunami in the Indian Ocean. 
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Figure 11. SLR and Crossover RMS for GFO POE arcs from January 4, 2000 to June 28, 2006. 


Table 7. SLR, Crossover, and Doppler RMS of fit for 
GFO POE Arcs 


Year 

Number 
of Arcs 

Average RMS of Fit 
SLR Crossovers Doppler 
(cm) (cm) (cm/s) 

2000 

67 

4.68 

8.41 

1.74 

2001 

60 

4.70 

8.64 

1.93 

2002 

66 

5.39 

8.12 

2.10 

2003 

63 

4.45 

7.12 

1.93 

2004 

62 

4.49 

6.80 

1.75 

2005 

71 

3.26 

6.57 

1.90 

2006 

34 

3.18 

6.58 

1.54 

All 

423 

4.37 

7.51 

1.89 


Table 8. Orbit Overlap Statistics for POE arcs from January 4. 
2000 to June 28, 2006 


POE Orbit Statistics (cm)* 
Mean Median Minimum 

Maximum 

Radial 

1.84 

1.39 

0.14 

15.39 

Cross-track 

11.56 

6.15 

0.57 

569.56 

Along-track 

13.50 

8.48 

0.83 

136.82 


* The orbit statistics are computed over 316 overlapping POE arcs. 
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Figure 12. GFO orbit error relative to TOPEX/Poseidon inferred from the mean RMS of the GFO sea surface variability: 
Top figure, before empirical correction; Bottom figure after empirical correction using TOPEX 


The precision orbit lies at the heart of the altimeter measurement, providing the reference frame from 
which the radar altimeter measurements are made. The NASA GSFC Planetary Geodesy Laboratory pro- 
duces a daily GFO orbit with a latency of less than 24 hrs and a precise orbit with a latency of three to 
four weeks. The modelling for the MOE and the POE orbits incorporates improvements to the gravity field, 
the nonconservative force model, and to the measurement model. After many improvements, the expected 
radial accuracy is approximately 5 cm for the POE, and between 10 to 20 cm for the daily MOE. Further 
improvements are possible for the MOE and POE orbits, for example using the GRACE generation gravity 
models, improved ionosphere corrections from the GPS Ionosphere Models, and more sophisticated handling 
of the nonconservative forces. 

An important lesson from GFO is that altimeter missions should always carry multiple means of tracking. 
Not only can they serve as a backup if one system fails, but the multiplicity of data types allows us to directly 
intercompare orbits computed using different data and verify the orbit quality, as we have done so successfully 
on TOPEX/Poseidon and on Jason-1. 4,5 
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